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Abstract 

We propose a method for reconstruction of the optical potential from scattering data. The algo- 
rithm is a two-step procedure. In the first step the real part of the potential is determined analytically 
1 via solution of the Marchenko equation. At this point we use a diagonal Pade approximant of the 

corresponding unitary 5-matrix. In the second step the imaginary part of the potential is determined 
via the phase equation of the variable phase approach. We assume that the real and the imaginary 
parts of the optical potential are proportional. We use the phase equation to calculate the propor- 
tionality coefficient. A numerical algorithm is developed for a single and for coupled partial waves. 
The developed procedure is applied to analysis of 1 S NN, 3 5£>i NN, P31 ir~N and 501 K+N 
data. 
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for the solution of the inverse-scattering problem at fixed angular momentum exist, and we have several 
of excellent reviews on the subject The direct application of these approaches to construction of 

local two-body potentials which are phase equivalent to the effective potentials occurring in theories 
describing reactions of composite particles is impossible. For, such potential must be complex in order 
to reproduce the loss of flux above the inelastic threshold. But it must reproduce the real phase shifts 
below this threshold and must be real itself. These requirements are incompatible for potentials being 
energy independent by construction. For very low threshold these approaches are applicable and produce 
energy independent complex potentials [oj Hy . In the general case the empirical energy dependent optical 
potentials are usually inferred by fitting of the parameters of an assumed analytic potential [7| |S] . This 
approach has two major shortcomings: a complexity and inconvenience of fitting simultaneously many 
nonlinear parameters; and lack of correlation of the parameters obtained at various energies. 

In this paper we develop an inversion method that is free of these shortcomings. The method is based 
on a fixed-Z inverse scattering theory and on a special parameterization of the optical potential. The 
proposed procedure is a two-step process. In the first step the phase shift data are used to determine a 
real potential via solution of the Marchenko equation. At this point we use a diagonal Pade approximant 
[M/M] of the corresponding unitary 5-matrix. In the second step the imaginary part of the potential is 
determined via the phase equation of the variable phase approach |13j . We assume that the real and the 
imaginary parts of the optical potential are proportional. The value of the proportionality coefficient is 
predicted by the phase equation and is refined by the iterative algorithm. We develop this method for 
single and for coupled partial waves. The whole procedure is applied to analyze 1 So NN, 3 SD\ NN 
data and to analyze P31 n — N and 501 K + N data. These analyses demonstrate that prediction for 
the proportionality coefficient from the phase equation is very close to a precise value that reproduce the 
experimental loss of flux. 

The plan of the paper is as follows. In Sect. 2 we describe the inverse scattering techniques based on 
the Marchenko integral equation. The used diagonal Pade approximants of the corresponding 5-matrix 
allow an analytical solution of the Marchenko integral equation [21 0] . For single partial wave the general 
solution was presented in [SJ|S]. We present a solution for coupled partial waves. These techniques 
produce real local potentials from phase shift analysis data. In Sect. 3 we consider the phase equation. 
We investigate how the 5-matrix is changing with certain change of the potential. This consideration 
shows advantages of proportionality of the real and imaginary parts of the optical potential. In Sect. 4 
the feasibility of the method is shown in the examples of analyses of NN, ir~N and K + N scattering 
data. 
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2 Inversion algorithm 



The Marchenko inverse scattering theory is viewed in detail in Refs. I2HH1E]- We shall, therefore, only 
briefly describe this formalism. 

The input data of the Marchenko inversion are 

{S(q),(0<q< +00), q~j, Mj, j = l,...,n b ], (1) 

where S (q) is the scattering matrix dependant on the relative momentum q, q 2 = Em, eft = mEj < 
0, Ej is the energy of the j-th bound state, so that iqj > 0, m is the reduced mass. The Mj matrices 
give the asymptotic behavior of the corresponding normalized bound states. 
We proceed from the Marchenko equation for single channel 

+00 

F(x,y) + L{x,y) + J L(x,t)F(t,y)dt = Q, (2) 

X 

where the input kernel is given by 

1 T 

F(x,y) = — / h+ (qx) (I - S (q))h+ (qy) dq + J^ Mfh+ (iq jX ) h+ {iq.y) , (3) 

00 i=i 

fif (z) are the Riccati-Hankel functions. 

The output kernel L(x,y) gives the reconstructed potential 

V W = -^. (4) 

This local energy independent operator V (r) links the Marchenko equation (0 and the radial Schrodingcr 
equation of a fixed angular momentum, 



d 2 1(1 + 1) T T { \ 

dr z r z 



^(r,q) = q 2 ip(r,q). (5) 



The scattering matrix S (q) , matrices Mj and energies Ej are the output data of the direct scattering 
problem associated with the Schrodinger equation (J5J ■ 

It has been known for several decades now that S matrices rational in q (ratio of polynomials) 
correspond to potentials known as Bargmann potentials expressible in terms of the elementary functions 
[3 El El- Such fraction may have the same truncated Taylor series as the S matrix it represents. It 
is then called a Pade approximant. Conjectures and theorems concerning the convergence and analytic 
continuation properties of Pade approximants are collected in |1(J) . For single partial wave the general 
solution of the Marchenko equation via Pade approximant of the S'-matrix was presented in jS] and in 

We shall, therefore, only present it briefly and turn to the case of coupled partial waves. 

A diagonal Pade approximant of the S'-matrix is given by 

s{q)=e , s= mzim (6) 

72 (q) + ifi (<?) 

fi (q) and f% (q) are an odd and even polynomials of q, which do not turn to zero at the real axis 
simultaneously. 

This approximant leads to the following expression for the phase shifts 5 (q) 

*<-«(.»-&$ m 

Inasmuch as lim^oo 5(q) ~ 1/q for regular potentials, it is evident that degree of fi(q) must be less 
than degree of /2(g) by 1. Let us select N discrete momenta % such that the corresponding 8 (%) = Si are 
known. Use of these values in eq. transforms the latter into a set of inhomogeneous linear equations 
from which N coefficients of the polynomials fi (q) and fa (q) can be determined. This is a usual strategy 
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[5], but since any set of qi is experimentally limited from above (qi < q ma x) there is some uncertainty 
in determination of S(q). Even if a set of qi is dense and the agrement between the data and the used 
approximant is excellent the arbitrary behavior of S(q) above q m ax guarantees that a solution of the 
inverse problem is arbitrary as well. We assume that about and above q ma x the true 5(q) depends only 
slightly on details of the potential and does not depend on its asymptotic at r > r max . Then we take a 
model potential (in our calculations V mo dei(r) = Aexp(— br)) and fit parameters (A and b) so that 

^modelis^ qmax) ~ qmax) > (8) 

here signs « mean that about q = q max the chosen model potential gives a phase curve which goes inside 
error bars. This means that we take 5 mo dei{q) as an asymptotic for S(q) when q > q m ax- In the segment 
[0, q m ax] the 5(q) is defined by some spline that approximates the data points {qi, Si}. In this way we may 
control the Pade fit in the line segment [0, Qmax], where Qmax is arbitrary and is fixed by convergence 
of the whole inversion procedure. The needed accuracy of approximant is attained by increasing of N 
which in turn defines degrees of fi(q) and f2(q)- 

Approximant J2J leads to a degenerate input kernel F (x, y). We calculate the integral in eq. JSJ using 
the residue theorem. For approximant © the result of the integration is 



+ £ Mfh+ (q t x) h+ %y) 



F (x, y) = % £ Res [h+ (qx) (I-S (q))h+ (qy)] 

i=l 

"pos Kb 

= £ hhf {fiix) h+ (fry) + £ Mfhf (q t x) h+ (q iy ) = £ b,h+ {ftp) h+ (fry), 

2—1 i—1 j — 1 



(9) 



where /3, (i = 1, ...,n poa ) are all S'-matrix poles with > 0, /3 = {(3\, ..fi n , <fi, ■■•3n i ,}, n = n pos + n&. 
We assume that all poles are of first order so that 



Res [h+ (qx) (I-S (q))h+ (qy)] = 2iRes [h+ (qx) M ^Sl_h+ (qy) 



<J=ft (10) 



here we have denoted f[(q) — dfi(q)/dq, (i = 1, 2). 

In this case the input kernel of eq. is a degenerate one as well as its output kernel 

n 

L(x,y) = ^Pi(x)h+((3 iy ), (11) 

i=l 

where Pi (x) are unknown coefficients. Substitution of 10 and into |J2} yields 

/ +?° \ 

£ h+ (fry) hh+ (J3 t x) + Pi(x) + h Pk (x) J h+ (J3 k t) ft+ ([3 z t)dt = 0. (12) 



i=l \ fc=l 



Linear independence of the hi (fry) implies that 



b l hi(l3 t x) + P l (x) + b l J2P k (x) j hi(/3 k t)hi(pit)dt = 0, (13) 

X 

or 

n 

J2^ik(x)P k (x)=D i (x) (i = l,..,n), (14) 

= -bihi (frx) a 

we have 



k=l 

H ' " 



where Di (x) — —bihi (frx) and after applying Riccati-Hankel integration formulas (see Appendix) in 



\ + hx((hi (frx)) 2 -hi^(f3 2 x)hl fl (frx))/2, for i = k 
-bi — — ■ — , for i k. 
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The functional coefficients Pi (x) are defined by I jl 4|) 

P l (x) = {A- 1 D) i . (16) 

Finally we derive L (x, y) and the potential V (r) from Ijll ll and 

In case of two coupled channels we present only sketchy derivations because of their awkwardness. In 
this case the system of the partial Schrodinger equations is 

^> =(*<;> Yrtt ) , <i8) 



V T (r) V 2 (r) 

where V\ (r), V2 (r) are potentials in channels 1 and 2, Vr (r) is potential coupling them, Xi( r ) ancl X2(^) 
are channel wave functions. 

By analogy with © we approximate the 5-matrix by the following expression 



qf \ _ f exp (2i5i) cos2e i exp (t ((Si + <5 2 )) sin2e \ _ 

I zexp (t ((Si + ^2)) sin2e exp (2162) cos2e y 



/« 1) ( 9 )+ l /< 1) ( 9 )y' (/< 12) ( 9 ))V(/< 12) ( 9 ))^ (/< 12) ( 9 ))%(/ 1 (12) (?))' .7=1,2 / a U, ( 9 )+.A u) ( 9 ) 

V (/i 12 ^))^ 12 '^)) 2 iii )2 / a 0) (9)-HA 0, («) U a) W+»A (2 W (/< 12) ( 9 )) 2 +(/ 1 (12) ( 9 ))^ J 

(19) 

This is again the most general Pade approximant for the S'-matrix. It was used in jl 1| in an other 
form, but corresponding analytical solution of the inverse scattering problem was not presented. 
The coefficients of this Pade approximant are determined from the equations analogous to l|20|) 

t j.m) = fpn, j=1 , 2 (a, 

and 

*ff (£(<?)) = TTY- (21) 

/2 (?) 

The generalized Marchenko equation for coupled channels formally has the former view |12j 

+00 



L(x,y) + F(x,y) + J L (x, t) F (t, y) dt = 0, (22) 

X 

but functions involved are matrices (2 x 2) 

+ °° n b 

F(x,y) = ±- I H(qx)[I-S(q)}H(qy)dq + J^H(p i x)M i H(l3 i y), (23) 

where 

»«-(£%(.,)•'-(;:)■ <«> 

Insertion of (|I9|I into (|23[) and applying of the residue theorem yields 

n pos n h 

F{x,y) = iY. Res [H (qx) (I - S (q)) H (qy)]\ q=/3i + £ H [fax) M?H (/%) = 



i=l 

,(2) 



(25) 



Y,H(p i x)Q\H{p i y)+ J2 xH>(p i x)Q 2 i H((3 i y)+ £ H ((3 t x) Q^H' ((3 t y)y, 

i—1 i—1 i—1 
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where /3j (i — 1, ...,n poa ) are all S'-matrix poles with > 0, /3j (i = 1, ...,rt^) are poles of the second 
order, /3 = {/3i, ~P n v) , ■■Pn po .> QU >»Qn t }, n = n pos + n h , 



it i ( \ _ ( dhf (x) jdx 

(X) ~ \ 1 dh+ (x)/dx 



We note that there are poles of the first as well as of the second order in the diagonal matrix elements 
and there are poles only of the first order in the off-diagonal matrix elements. Poles of the second order 
in the diagonal elements are poles of the first order in the off-diagonal matrix elements, and they must 
be enumerated twice. Q\ (j = 1, 2) are constant matrices they are trivial but cumbersome therefore, we 
do not give them. 

We solve eq. <|22[1 using substitution 

n n 

l (x, y) = J2 p i W H (ftw) +J2 N i ^y H ' - ( 26 ) 

where Pi (x), Ni (x) are unknown functional (2 x 2) matrix-coefficients. Linear independence of the 
H (J3%V) an d yH' (P%y) implies that 



E p i 0*) Q% (*) +T, N i {x)Q% (x) = H (p jX ) Q) + xH> (p jX ) Q) 
E k (z) Q& +E p i (*) = H (fts) Q 2 



(27) 



where 



(a;) = + / H (^t) H %t) dt x Q) + j tH H' ((ijt) dt x Q] 



+ OC 

[x)= J H [fat) H (pjt) dt x Q 2 (28) 



+oo +oo 

%{x)= [ tH' (Pit) H (fijt) dtxQ)+ ( t 2 H' (f3it)H f (fytydtxQ? 



X X 



Q?- (a:) = ISij + J tH' (^t) H (Pjt) dtxQ'j, 

X 

Integrals of expressions i|28|) are presented in Appendix. Matrix equations (|27H can be trivially reduced 
to scalar linear equations. Having solved this linear equation system we get the sought-for potential from 
and 0). 

The multichannel generalization is trivial. 



3 The optical potential 

In this section we consider changes of S'-matrix that are induced by certain transformation of real 
potential. 

First we consider the one channel problem. 

The phase equation for the initial potential V° (r) obtained by some inversion procedure (from 
Marchenko equation in our calculations) is 

oo 

^ (0) =-- J V {0) (r) D? (qr) sin 2 (S t (qr) + <S (0) (r)) dr, (29) 
o 
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where £>i (z) and Si (z) are Riccati-Bessel amplitude and phase correspondingly ^3] 



D l (x)=^jf(x)+nf(x), (30) 

Si(x) = — axctan(ji(x)/ni(x)) (31) 
Let us consider the complex-valued potential V^ 1 ) (r) obtained from V° (r) by transformation 

V {1) (r) = (1 + to) V {0) (r) , (32) 

where a is some real parameter. Such parametrization was used in |14| but without analysis (a was 
fitted). Evidently the phase equation for this potential is 



<5 (1) = -- (1 + iot) { V<°> (r) D? (qr) sin 2 ^(qr) + (r)) dr. (33) 



q 

o 



From eqs. and lO we get 

tfW - (l + ia)ff<°> = 

/V<°> (r) £>f (gr) (sin 2 pj(gr) + <5« (r)) - sin 2 (§i(qr) + 5^ (r) 



, 1+ia 
9 





1+ig 
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dr = 



(34) 



/V<°> (r) A 2 (<7 r ) sin (?k{qr) + 5^ (r) + 5^ (r)) sin (5^ (r) - <y(°) (r)) dr 



For smooth enough potentials the right side of eq. I|34|l rapidly decreases comparing with S^ ' and 
because there is a rapidly oscillating around zero function under the integral in (|34|) (underlined). 
Its frequency behaves as 2q for big q (see l^lf) ). Then as the first approximation we may take 

S {1) « (1-Ha)<5 (0) =S R + iS I . (35) 

For inelastic scattering the S'-matrix is expressed through the real inelasticity parameter p and the 
real phase shift 5 

S = cos 2 (p) e 2lS = e 2 ^ 5R+lS '\ (36) 

so we easily arrive at 

5 R = 5^5^\ (37) 

cos 2 (p)«e- 2a5<0) , (38) 

whence it follows that a8 > 0. The formula l|38[l allows to calculate the parameter a from the known 
values p and 5^ ~ S. 

Eqs. I|34l35jl imply that 

DO OO 

F (0) (r) Df (qr) sin 2 (l,(c?r) + <5 (0) (r)) dr ~ / (r) £ 2 (gr) sin 2 (<S ; (gr) + <5« (r)) dr, (39) 
o o 
Consideration of the coupled partial waves is more complicated. The initial real potential is 

/ T/ (0) T/ (0) \ 

^w-(^^J- (40) 

The equations for eigenphases and mixing parameter of potential (|40|l are J3| 

(41) 
(42) 
(43) 





_ r(0) 
— J ll ' 




, j(0) 
^ i 13 


,(1) 

°2 


_ r(0) 
— J 21 ' 




, r(0) 
^ -'23 




- / (0) ^ 

— ^31 " 


-4^ 


AO) 

' J 33 ) 
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where 



r(0) 



- f drV{ 0) {r) cos 2 eW(r)Dl(qr) sin 2 (6 h {qr) +Sf\r)) 



t[°) = -- [ drV}°\r) sm 2 eW(r)Dl(qr) sin 2 (d l2 (qr) +6{ 0) (r)) 



T (P) 

'13 



- [ drV^\r)sm2e( \r)D h {qr)sm(5 h (qr) + 5 i °\r))D ll {qr) (44) 



4i ] = -- / drVf 0) (r) sin 2 e(°\r)Dl(qr) sin 2 (5 h (qr) + S^°\r)) 



I<°> = -1 /"dryi 0) (r)cos 2 e( )(r)A 2 2 (gr)sin 2 (4(?r)+4 0) (r)) 



^ = -- / dr^ 0) (r) sin 2e(°)(r)A 2 ( (Z r)sm(4( (? r) + 4° } (r))/)^ (gr) sin^ (qr) +^ 0) (r)) (45) 



r(0) _ 
J 31 — 



?(0) 2 ,1 W ^(O) , » V i 0) (gO sin (4 (g*-) + ^ (?)) Mk (gr) + 4° } (0) 



2<?7 sin(^(r)-4 0) W) 



r(0) 
'31 



sin2e(°)(r)dr 
2? 7 sin(^ 0) (r)-4 0) (r)) 



F 2 (0) (r)A 2 2 (<zr) sin^gr) + 6\ v '(r)) sm(S h (qr) + 8\ u >(r)) 



(0), 



r(0) _ 1 fV^°\r)D h (qr)D h (qr)dr 



'31 



29 7 sin(<5i 0) (r)-4 0) (r)) 



2e(°)(r)sin(^ 1 (gr) + ^ 0) (r))sin(4(gr) +4 0) (r)) 



1 (cos 2e(°) (r) - l) sm(8 { ° } (r) - 5^ (r)) sin (\ (gr) - <5, 2 (qr) 



(46) 



By analogy with the one channel case the following generalization for the optical potential is derived 



v( i) {r)= I (l + ia{)V^ (l + ia 3 )V 5 



-(f)) 



v (l + ia 3 )V± 0) (l + ia 2 )V 2 {Q) J ' 
Evidently the phase equations for this potential is 

S[ l) = (1 + iai)I$ + (1 + ia 2 )I$ + (1 + ia 3 )l£> 

^23 
^33^ 



(1 + iai)/^ + (1 + ia 2 )I$ + (1 + 103)/. 
(1) = (1 + iai)4\ } + (1 + ia 2 )42 } + (1 + iQ! 3)/ ; 



(i) 



(47) 



(48) 
(49) 
(50) 



Integrals 7^ are defined as 7^ in I|44I46|I but through 5^(r), 6^ (r) and e^^(r) instead of <5^(r), 



r (o) • 



4 0) (r) ande(°'(r). 

Evidently we cannot consider I|48I5U|I in a manner like (|34fl . But we assume that 

^ = ^+ E ^ 

i,i=l,2,3 



(51) 
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where 

E °(^ ) )«^ ) ' far i,j = 1,2,3. (52) 

i,j=l,2,3 

This assumption can be considered as a generalization of 1)39(1 . It is hard to prove in the general case, 
but our calculations show that this is true at least in case of 3 SDi NN scattering. 

Eigenphases i — 1, 2 and mixing parameter e^°) are real and they define a unitary S'^-'-matrix 



/ cos 2 e(°)e 2 4 0) + sin 2 e(°) e 2 4 0) cose^) sine^ (e 2 <' 

cose(°) sine(°) (e 2 <' - e 2 ^"') sin 2 e(°) e 2 ^°' + cos 2 e^e 2 ^ 



e 2i5 2 



(o) 



V 



(53) 



Eigenphases <^ and mixing parameter e' 1 ) are complex but they define S^-matrix in the regular 



way 



/ cos 2e(i) e 2'«! 1) + sin 2 e( 1 )e 2 ^ 



cos e^ 1 -* sin e( 



1) LttS™ 



y coseW sine^ (e 21 ^ ' 



sin 2g(l) e 2*<5« +cos 2g(l) e 2»«( 1 



(54) 



The complex eigenphases <5 4 - and mixing parameter e^ 1 -* are defined by the experimental S'-matrix (5 
S^). A direct consequence of (|48I51|) is 



mP=Sf\ i = l,2; Sfe^W ). 



^E^ i = 1 = 2 > N x) = E«i4 0) - 



(55) 
(56) 



We may calculate coefficients Jy ~ if^ — 1,2,3) using 1(441461) . A simpler method is to use the 
following implication of 148151( 1 



i = l,2; I. 



3 ^ ~^~9a. 



(57) 



Calculated coefficients may be checked by ((41143(1 . Next, we calculate on (i = 1, 2, 3) from 1(48150(1 . 



4 The optical potentials 

We apply the developed method of inversion to analysis NN, tt~N and K + N data up to energies where 
relativistic effects are essential. We take into account these effects in the frames of relativistic quantum 
mechanics of systems with a fixed number of particles. The review of this approach can be found in |15) . 
Here we give only some extracts of it. 

The relativistic quantum mechanics of systems consisting of a fixed number of particles is based on 
the conjecture that the number of particles is constant at not very high energies and on the assumption 
that the group of invariance for the system under consideration is the Poincare group rather than the 
Galilei one. A system of two particles is described by the wave function, which is an eigenfunction of the 
mass operator. In this case we may represent this wave function as a product of the external and internal 
wave functions (161 117( . The internal wave function \ is also an eigenfunction of the mass operator and 
satisfies the following equation 



\Jq 2 + m 2 + ^q 2 + m 2 + V» 



X = Mx, 



(58) 



where Vint is some interaction operator acting only through internal variables (spins and relative momen- 
tum), q is a momentum operator of one of the particles in the center of masses frame. Rearrangement of 
j5B) gives 

[q 2 + 2mV] X = q\, (59) 



8 



where 



4 2 4M 2 

to is taken as a nonrelativistic reduced mass 



9 M 2 to? + to? (to? — to?,") 2 . s 

1 2 = -. ■ (60) 



_mima_ 

TO,i + TO2 



1/ is an operator acting like V^„t only through internal variables. 
In case of two particles with equal masses toi = mi = 2m 

9 2 = ^-4to 2 . (62) 

Eq. I|59l) is identical in form to the Schrodinger equation. The quasicoordinate representation cor- 
responds to the realization q = —i-§^, V = V(r). In [23 we showed that this formalism can be easily 
generalized for the case of inelastic channels, particularly it allows to take into account isobar channels 
in NN scattering. This formal coincidence allows us to apply our inversion algorithm. 

We applied the described algorithm of inversion to reconstruction of the nucleon-nucleon potential. 
As input data for this reconstruction we used modern phase shift analysis data up to 1100 MeV for 3 SZ?i 
state and up to 3 GeV for 1 So state of nucleon-nucleon system OI2J3]. The deuteron properties were 
taken from |21j . These data allow to construct nucleon-nucleon partial potentials sustaining forbidden 
bound states (Moscow potential introduced in (22). Parameters of forbidden bound states for the partial 
Moscow potentials were chosen to be equal to those of model potentials (see Sect. 2). In this way 
we constructed the NN optical potentials for 1 So and 3 SDi partial waves. These potentials describe 
the deuteron properties and the phase shift analysis data. The 1 5'o phase shifts of Moscow potential 
begin from tt. 3 So phase shifts of Moscow potential begin from 2ir. The mixing parameter ej of Moscow 
potential differs from that of traditional repulsive core potential by sign. The real parts of the constructed 
partial potentials are presented in fig. 1. The calculated values of deuteron properties are compared with 
the experimental data |21) in Table 1. Only three parameters are fixed as input data of inversion problem. 
These parameters are energy, A$ and r\dls- Figs. 1 and 3 demonstrate how changes of 5 influence the 
partial potential for 1 So wave. 

As another example of application we analyzed the modern P31 ir~N data up to 2 GeV and 5*01 
K + N data up to 1 GeV (20j and constructed the corresponding optical potentials. The real parts of the 
constructed partial potentials for P31 -k~N data up to 2 GeV and SOI K + N are presented in fig. 2. 

From eqs. I|38[) and I|56II57[I we calculated parameters a and cti, (i — 1, 2, 3) which define the imaginary 
parts of potentials. Our predictions were justified. Calculations with optical and real potentials (fig. 3,7,8) 
show the validity of l|37|) and JSJJ. The a's predicted by l|38|) may be improved by a simple numerical 
method. Predicted and improved values of a's are shown in fig. 4,5,6. In all figures "Calc. I" means 
calculations from predicted values of (|37|l . "Calc. II" means calculations from refined values. Fig. 5 
shows that calculation of on, (i = 1,2,3) from (|56l and l57() does not require refinement because it uses 
more precise values of iff than those implied by assumption that iff = iff . 

In eq. (|36|l we use parametrization of the S-matrix from |23| whereas parametrization of partial wave 
analysis is based on type- if scheme P3] . For uncoupled waves, the S'-matrix is given by 

c _ l-K,+iK r 

~ 1+Ki-lKr ' 

(63) 

where K r = tanS, K$ — tan 2 p. 

So we had to recalculate data of hato S'-matrix, then into parameters of [2H| to get input data of 
inverse problem. Our results are presented in parametrization of |19j . 

All potentials and inelasticity multipliers (a's) may be downloaded from cite www.physics.khstu.ru 
in numerical form. 

5 Conclusions 

Let us summarize the results presented in this work. In the first place we mention a presented analytical 
solution of the Marchenko equation for coupled partial waves in case of diagonal Pade approximant of the 
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corresponding S'-matrix. The inverse scattering scheme at fixed angular momentum is used to construct 
a local real energy independent potential as a first step of our inversion procedure for single and coupled 
waves. Furthermore, we consider what changes of S'-matrix are induced by certain transformation of real 
potential. We have found out that certain simple transformation may have a negligible effect on phase 
shift but introduce a controllable inelasticity. This transformation does not change the real part of the 
potential but adds an imaginary part. As a result we get an optical potential with energy independent 
real part and energy dependent imaginary part. We apply this scheme to NN, n~ N and K + N scattering 
successfully. 

6 Appendix 

We present only nontrivial integrals of l|28l) . They can be derived from the recursion relations for the 
Riccati-Hankel functions and from known integrals |26| . 



h(x,fa,fa,l) = J h+(fat)h+(fa)dt = 



J h'+(fat)h + (fat)tdt 



Where fa k = (3? - 0* 



-x((hf(fax)) 2 - ft+ ! (ftz) fy^ (ftz))/2, for i = 

, for k. 



p t h+_ l (!3 % x)h+ (I3 k x)-f3 k h+ ^fffcx) 
Pit, 



-(x(h+(fax)) 2 +h(x,fa,fa,l))/(2fa), for i = k 
§| (fah + (fax) h+_ t (fax) - fahf^ (fax) h+ (I3 k x)) + 
+ ^h+_ 1 (fax)h+(fax)- 
-f-{fah\ + (fax)h+_^fax) - fah[t x (fax)h+ (fax)) , for i^k. 



(64) 



(65) 



J h'+(fat)h'+(fat)t 2 dt = 



h [fax, I) /ft, for i = k 
S^f- {fahf (fax) h+_ t (fax) - fahf (fax) (fax)) + 
(h+ (fax) h+_ x (fax) + faxhj (fax) h'+ x (fax) - 
-faxh'+ (fax) h+_ 1 (fax)) + 

+f~ k (K + (Pk*) h ti (Pi x ) ~ K + (Ax) K-i (Pkx) - 
-faxti+ (fax) ti+ x (fax) + faxh'+ (fax) h'+ x (fax)) - 
+ f| (h+ (fax) h+_ t (fax) - faxh'+ (fax) hf_ t (fax) + 
+faxhf (fax) h'+ x (fax)) , for i^k. 

I 2 ( Zj l) = ( lz 2 /2 - 3z/2 - 9i/2 + 1/z) exp(z2z) 

h (z, 2) = (-iz 2 /2 + 7z/2 + 49i/2 - 24/z - 24i/x 2 + 12/x 3 ) exp(i2z) 

h (z, 3) = (iz 2 /2 - liz/2 - 169z/2 + 171/z + AhOi/x 2 - 765/a; 3 - SlQi/x^ + 405/x 5 ) exp(z2z) 

I 2 ( z , 4) = (~iz 2 /2 + 21z/2 + 44U/2 - 745/z - 3510i/x 2 + 11835/x 3 + 
+28560i/a; 4 - 47880/a; 5 - 5040(h/x 6 + 25200/x 7 ) exp(z2z) 



(66) 



(67) 
(68) 
(69) 

(70) 



References 

[1] B. M. Levitan, Generalized Translation Operators and Some of the Applications (Davey, New York, 
1964). 

[2] Z. S. Agranovich and V. A. Marchenko, The Inverse Problem of Scattering Theory (Gordon, New 
York, 1963). 

[3] R. G. Newton, Scattering Theory of Waves and Particles (McGraw-Hill, New York, 1982). 

[4] K. Chadan and P. C. Sabatier, Inverse Problem in Quantum Scattering Theory, 2nd ed. (Springer, 
New York, 1989). 



10 



[5] A. Papastylianos, S. A. Sofianos, H. FiedeldeyAlt, E. O. Alt, Phys. Rev. C42, 142 (1990). 

[6] E. O. Alt, L. L. Howell, M. Rauth, S. A. Sofianos, Phys. Rev. C48, 176 (1993). 

[7] P.V. Green, K.W. Kemper, P.L. Kerr, K. Mohajeri, E.G. Myers, D. Robson, K. Rusek, and I.J. 
Thompson, Phys. Rev C53, 2862 (1996). 

A.Funk, H.V. von Geramb, and K.A. Amos, Phys. Rev C64, 054003 (2002). 

[9] Th. Kirst, K. Amos, L. Berge, M. Coz, and H.V. von Geramb, Phys. Rev. C40, 912 (1989). 

[10] G. A. Baker, P. Graves-Morris. Pade Approximants. Massachusetts: Addison- Wesley Publishing Co. 
1981 

[11] J. Wiesner E.P. Zhidkov, V. Lelek et al, Physics of elementary particles and atomic nuclei, V.9. Part 
3, 710 (1978). 

[12] M. Blazek, Matem.-Fys. Casopis, 13, 147 (1963); Commun. Math. Phys., 3, 282. 

[13] F. Calogero, Variable phase approach to potential scattering, Academic Press, New York and London, 
1967. 

[14] V. G. Neudatchin, N. P. Yudin, Y. L. Dorodnykh, and I. T. Obukhovsky, Phys. Rev. C43, 2499 
(1991). 

[15] B.D. Keister and Polyzou, Adv. Nucl. Phys. 21 , 225 (1991). 
[16] F.M. Lev , |nucl-th/9403222. 



[17] N.A. Khokhlov, VA. Knyr, V.G. Neudatchin, Phys. Rev. C68, 054002 (2003). 
[18] V.A. Knyr, N.A. Khokhlov, Phys. of At. Nucl., V. 67, N5, 937 (2004). 

[19] RA. Arndt, I.I. Strakovsky and R.L. Workman, Phys. Rev. C62, 034005 (2000); R.A. Arndt, J.S. 
Hyslop and L.D. Roper, Phys. Rev. D35, 128 (1987); J Bystricky, C. Lechanoine-Leluc, F. Lehar, 
Eur. Phys. J. C4, 607 (1998); J. Ball, J. Bystricky, J.-M. Fontaine et al, Nuovo Cimento Alll, 13 
(1988). 



[20] R.A. Arndt, W.J. Briscoe, R.L. Workman, I.I. Strakovsky , http://lux2.phys.va.gwu.edu/ 

[21] R. Machleidt. Adv. In Nucl. Phys. V. 19, 189 (1989). 

[22] V.G. Neudatchin, I.T. Obukhovsky, and Yu. F. Smirnov, Phys. Lett. 43B, 12 (1973); V.G. Neu- 
datchin, I.T. Obukhovsky, V.I. Kukulin, and N.F. Golovanova, Phys. Rev. Cll, 128 (1975). 

[23] L. Ray, Phys. Rev. C35, 1072 (1987). 

[24] R.A. Arndt, L.D. Roper, R. A. Bryan et al, Phys. Rev. C28, 97 (1983). 

[25] R. A. Bryan, Phys. Rev. C39, 783 (1989). 

[26] I.S. Gradstein, I.M. Ryzhik, Table of Integrals, Series and Products, (Academic, New York 1980). 



11 



Table 1. The deuteron properties 





Exp. a 


Calculation 
with Moscow 
potential 


Energy (MeV) 


2,22458900(22) 


2,2246 a 


Q (Fm a ) 


0,2859(3) 


0,277 c 


A s (Fm- 1 /^) 


0,8802(20) 


0,8802 


r d (Fm) 


1,9627(38) 


1,956 


Vd/s 


0,02714 


0,02714 




0,857406(1) 


0,859 c 



a relativistic correction included; b Data are from [21] ; c Meson exchange currents are not included. 




Figure 1: Real parts of NN potentials. Solid line 1 S'o -one channel (inversion from new fit of fig. 2). 
Dash-three points line 1 5'o -one channel (inversion from old fit of fig. 3). Two bound channels long dashed 
line Vcs(r) ( 3 <Si), short dashed line Vcd(t) ( 3 £>i), dotted line V te ns{r)- 
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Figure 3: Phase shifts and mixing parameter. Results for optical potentials are indistinguishable from that 
as for real potentials. Solid lines, reconstructed Moscow potential; dashed lines reconstructed repulsive 
core potential. The phase shift analysis data are from ^JlJOj- For S waves the original data set from 
[191 [SO] is raised 180 degrees up. To leave the 5-matrix unchanged we must then change the sign of the 
mixing parameter ei for the Moscow potential. SP00 and SP04 data are from |2(J| . 




Figure 4: Left: The inelasticity parameters. The phase shift analysis data (circles) are from ^5)- Right: 
parameter a. 
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Figure 6: Parameters cti for 3 SDi NN potential predicted by eqs. ljSl)|) and l(57|) . Solid line, a\] dashed 
line with circles, a^; dotted line, a^. 
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E |ab , MeV 



Figure 8: T-matrix and a for SOI K + N, data are from |2U) . 
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